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Is cosmography a useful tool for testing cosmology? 
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Model-independent methods in cosmology has become an essential tool in order to deal with an 
increasing number of theoretical alternatives for explaining the late-time acceleration of the Uni¬ 
verse. In principle, this provides a way of testing the Cosmological Concordance (or ACDM) model 
under different assumptions and ruling out whole classes of competing theories. One such model- 
independent method is the so-called cosmographic approach, which relies only on the homogeneity 
and isotropy of the Universe on large scales. We show that this method suffers from many short¬ 
comings, providing biased results depending on the auxiliary variable used in the series expansion 
and is unable to rule out models or adequately reconstruct theories with liigher-order derivatives 
in either the gravitational or matter sector. Consequently, in its present form, this method seems 
unable to provide reliable or useful results for cosmological applications. 

PACS numbers: 98.80. Jk, 04.20.0, 04.50.Kd, 95.36.+X 


I. INTRODUCTION 

The phenomenon driving the acceleration of the Uni¬ 
verse remains one of the biggest puzzles in fundamental 
physics today. In order to deal with this problem many 
models have been proposed, ranging from assuming an ef¬ 
fective perfect fluid with negative pressure, dubbed dark 
energy, modifications of General Relativity (GR) on large 
scales or even breaking the Copernican Principle (see [1] 
for reviews). As the number of possibilities increases, a 
useful tool is to use frameworks, which are able to en¬ 
compass a large class of models or theories and where 
observational scrutiny is able to rule out whole classes 
of possibilities without having to analyse each one on 
a case-by-case basis. Such model-independent methods 
have been applied to infer the dark energy equation of 
state and to reconstruct classes of dark energy theories. 
Among them, the cosmographic approach [2] , which relies 
on the Copernican Principle (leading to the Friedmann- 
Lemaitre-Robertson-Walker (FLRW) metric), and the 
expression of the scale factor as a function of an aux¬ 
iliary variable - either time or redshift for instance - have 
been used to look for deviations in the standard ACDM 
model or to reconstruct models for dark energy and the 
action for modified gravity [3-8]. 

In this paper, we discuss several limitations of cos¬ 
mographic analyses to constrain various models as well 
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as their ability to reconstruct theories from cosmological 
data. 


II. COSMOGRAPHY APPROACH 

The cosmography approach starts by defining the cos¬ 
mographic functions, 

a a . aW 

H ~ a ’ Q ~ aH 2 ’ J — aH 3 ’ S ~ aH A ’ ~ aH 5 ’ 

(!) 

where a is the scale factor, the dots represent cosmic time 
derivatives and a^ holds for the n th time derivative of a. 
Evaluated today, those functions define the cosmographic 
parameters Hq, qo, jo, so and lo, which are the Hubble 
constant, the deceleration parameter, the jerk parameter, 
the snap parameter and the lerk parameter, respectively. 
Then, by using the definition of the redshift 1 + z = -, 
the Hubble parameter H can be expanded in powers of 
the redshift around z = 0 as H = Hq+H^z+^y^z 2 +..., 
which converges for |;z| < 1 and where the Hubble param¬ 
eter derivatives, denoted by a z subscript, are evaluated 
today (subscript 0). They can be expressed in terms of 
the cosmographic parameters as follows (throughout this 
paper, let us assume that the curvature density parame¬ 
ter Dfc is negligible): 

H z o/Ho = 1 + qo , H zz o/Ho = —q% + jo , 
Hozo/Hq = 3(?q(1 + go) ~ jo(3 + 4g 0 ) — So , 

Hi z o/Ho = —3gg(4 + 8go + 5go) + jo (12 + 32go 

+25go — 4jo) + sq( 8 + 7go) + lo ■ (2) 
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In principle the parameters (2) can be fitted with obser¬ 
vational data, leading to model-independent constraints. 
As an alternative to the independent variable z, the above 
expansion may be expressed in terms of a more conve¬ 
nient variable that ensures the convergence of the series 
for the whole past history of the universe [4]: 


In order to illustrate the differences between the z and y 
redshift measures, let us use mock data generated from a 
fiducial spatially flat ACDM model, assuming the same 
redshift distribution as the Union 2.1 catalogue [9], with 
errors of magnitude a^ = 0.15. In this paper we restrict 
ourselves to the y and z parameters and supernova type 
data. Although many different observables have been 
used to constrain cosmographic parameters, such as H ( z ) 
data, baryon acoustic oscillations, gamma-ray bursts, an¬ 
gular distances to galaxy clusters [10], it is easier to un¬ 
derstand what is going on with only one observable. Thus 
we leave combined probes analyses for future work. 

Table I: Coverage test for 0i and 02. Refer to the bulk of the 
text for further details. 
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By generating 100 simulations based on this model, 
with = 0.3 and Hq = 73.8 km s _1 Mpc -1 , we 
then constrain the cosmographic parameters (1) by using 
both the y and z variables and then derive the poste¬ 
rior probability of the parameters. We have considered 
two sets of parameters: = {H 0l q 0 , j 0 , Sq} and 0 2 = 

{Hq, go, jo, Sq,Iq}, where H 0 is analytically marginalized, 
since supernovae do not provide constraints over it. We 
have sampled from the posterior and obtained the con¬ 
straints through a Monte Carlo Markov Chain (MCMC) 
using emcee [11]. The constraints coming from the mock 
data sets allows us to infer how frequent the true values 
are within the 1, 2, and 3cr confidence regions, thereby 
providing a way to test whether there was any bias for a 
parameter for each of the two expansions performed (y 
or z). 

Table I shows the number of times the true param¬ 
eters were inside the confidence region bounds for the 
first set of parameters in the y and z parameterisations. 
We would expect the true value to lie within la, 68% 
of the times and 2cr, 95% of the times. As shown, the 
z-parametrisation gives well-behaved coverage results. 
Only so overestimates the errors. On the other hand, 
the y-parametrisation gives completely biased estimators, 


where only a few simulations are within la and most 
of them lie at 3 or more a. Interestingly enough, this 
happens despite the fact that the y-parametrisation has 
much larger errors than the z-parametrisation. 

In the Figure 1 we have included a comparison between 
the exact ACDM Hubble parameter evolution and its cor¬ 
responding expansion in terms of the cosmographic pa¬ 
rameters following the 0i and 02 sets. For the z-redshift 
analysis, at z ~ 1 the differences are around 0.6% for 
both orders of the expansion whereas for the y-redshift 
analysis at a similar era y ~ 0.5 (z ~ 1) the differences 
yield around 9% while considering four parameters and 
even larger for three cosmographic parameters, showing 
the unexpected inadequacy of the y-redshift in a cosmo¬ 
graphic analysis. Figure 2 then depicts a given realisation 
of our simulations, so the comparison of the size of the 
error can be made. The trends are as follows: the y- 
parametrisation provides much bigger errors, biasing go 
to smaller values and jo and sq to greater values than the 
true ones. 

In addition, while considering higher order terms in the 
expansion (see Table I), the inclusion of a new parameter 
leads to an extreme overestimate of the errors for both 
variables, i.e., the values within la for almost all realisa¬ 
tions, where only go is less affected. These results show 
that, despite y might be preferable from a theoretical 
point of view, it is clearly not suitable to derive cosmolog¬ 
ical constraints, since the fittings are completely biased. 
One additional avenue that could be explored would be 
to derive a pivot redshift, around which the expansion is 
optimal. In theory, this could shed light on whether or 
not the behaviour we observe remains. Clearly, such a 
study depends strongly upon the redshift distribution of 
the sample under consideration and it is not within the 
scope of this paper. More to the point, this kind of study 
would be very helpful in order to understand how the 
redshift distribution of the sample influences the results, 
which might be important for some ongoing and future 
surveys as DES and LSST. 


III. TESTING ACDM 

One of the predictions of the flat ACDM model re¬ 
garding cosmography is that jo = 1. Consequently, we 
can easily check the consistency of the model by analyz¬ 
ing the posterior of this parameter. However, how the 
value of jo changes with the assumptions of the cosmo¬ 
graphic analysis (for instance the order of the expansion) 
and how the constraints can be compared to the ones 
obtained when fitting the free parameters of the model, 
directly determine the usefulness of the test. 

In order to address those questions we generated some 
mock realizations of data for a flat XCDM model, i.e., 
a model where the dark energy equation of state is 
given by w = px/ px, with H m = 0.3 and w = —1.3, 
which gives jo = 1.945. Fig. 3 (left and middle) dis¬ 
plays the constraints of a typical realization consider- 
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Figure 1: Left panel: Hubble expansion rate for exact ACDM model with = 0.3 and its corresponding cosmographic 
expansion with the y-parametrization with three {go, Jo, so} and four {go, Jo,so,fo} cosmographic parameters. Right panel: 
The same for the z-parametrization. 




Figure 2: Left panel: Constraints to the first set of cosmographic parameters with the j/-parametrization. Right panel: The 
same for the z-parametrization. 


ing the cosmographic expansion 0 \ (fourth order), Q 2 
(fifth order) and the exact XCDM model, H 2 /Hq = 
D m (l + z ) 3 + (1 — f2 m )(l + z) 3 ( 1+w \ As shown in Fig. 3, 
fitting the XCDM model directly provides smaller errors 
than the cosmographic approach, and it can spot devi¬ 
ations from ACDM with less effort. Moreover, the con¬ 
straints on w show a clear signature of w ^ — 1. On 
the other hand, the order of the cosmographic expansion 
decisively affects the posterior constraints on the cosmo¬ 
graphic parameters: there is some evidence of Jo 7 ^ 1 
when considering 6 1 but such evidence disappears when 
62 is assumed. 


IV. CONSTRAINING DARK ENERGY 
MODELS WITH COSMOGRAPHY ANALYSIS 


Recently it has been shown that cosmography can be 
used to reconstruct particular models for dark energy [ 6 ] . 
However, as demonstrated below, there are issues with 
this approach when it is applied to theories that con¬ 
tain higher derivatives, such as fourth-order gravity or 
Galileons. This is because the expansion of such mod¬ 
els give extra free parameters apart from the usual cos¬ 
mographic ones. We now compare several dark energy 
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models beginning with one without higher derivatives: 


\r- V{4>) +JZ m , 

. (4) 

where C m is the matter Lagrangian density, w(cj) is the 
factor that renormalises the scalar field </ and V(<f >) its 
potential. Note that for the sake of simplicity in the fol¬ 
lowing we use units where 87 tG = 1. Thus, the derivatives 
of the potential evaluated today, i.e., at redshift zero, can 
be written as 



Vo 

HZ 


v z0 

HZ 

V320 

“HT 


„ 3f l m 

= 2 -®- — , 

= 4 + 3®-i„-^ , 

= 4 + 8qo + jo(4 + 90 ) + so — , 

= Jo -lo- qojo(7 + 3go) - so(7 + 3q 0 ) 


(5) 


9Q m ■ 


where field equations in a spatially-flat FLRW were used. 

Assuming the universe today is well approximated by 
a ACDM model, we can write « 2/3(1 + q 0 ). Thus, 
in this case there is a one-to-one correspondence between 
the derivatives of the potential and the cosmographic pa¬ 
rameters. In this way, it is interesting to see how the 
method performs compared to other model-independent 
approaches. In Fig. 4 (left), the evolution of the scalar 
potential V((f>) and the Hubble parameter (right) are 
shown for the ler and 2 er regions stopping at third and 
fourth orders for a generic realisation of a spatially flat 
ACDM model. This can be compared with a Gaussian 
process regression (see Fig. la in Ref. [12]), where the 
errors are much smaller and no assumption regarding the 
model behaviour today was made. As can be seen, after 
z ~ 0.5 no useful constraints are derived for five pa¬ 
rameters, while for four parameters good constraints are 
obtained up to z ~ 1 . 

As an example, where higher-order derivatives are in¬ 
volved, let us now consider the case of f{R) gravity: 



+ £m 


( 6 ) 


The f(R) function and its derivatives with respect to the 
redshift z evaluated today lead to 


fo 

e m 

fzO 

6 Hi 

f2z0 


= -aq 0 + + 6/3 (2 + q 0 - j 0 ) , 


= a (2 + g 0 - jo) , 


(7) 


~~ 6/1 (2 + qo — jo) + 0 . [2 + 4go + (2 + qo)jo + so] 


In this case, there are two extra free parameters in com¬ 
parison with quintessence, namely 


d/ 
d R 


d 2 / 


R—Rq 


AR 2 


= ( 8 ) 

R—Ro n 0 


Previous works in the literature [7] fixed the values of 
a = 1 and fi = 0 a priori, such that the model coincides 
with General Relativity at z = 0. However, whenever 
d 2 //di? 2 = 0 either a singularity or instability occurs 
[13]. Apart from this theoretical shortcoming, it turns 
out that cosmological values a / 1 and fi ^ 0 may still 
produce viable cosmological models. Therefore, the naive 
assumption about the a and f3 parameters must be aban¬ 
doned and these two parameters should consequently en¬ 
ter in the analysis as free parameters. It follows therefore, 
that it is not possible anymore to have one-to-one cor¬ 
respondence between the /(A)-derivatives and the cos¬ 
mographic parameters. This means one should constrain 
the cosmographic parameters and use eqs. (7) and (8) 
in order to reconstruct / and its derivatives. As one can 
see, the data does not provide any constraints over a and 
/?, so either priors over these parameters or complemen¬ 
tary tests are necessary [14] . In order to illustrate the 
difficulties in defining sensible priors over the aforemen¬ 
tioned parameters, we generated mock data for the fol¬ 
lowing toy-model: f(R) = R + aR 2 + bR 3 with a = 2.81 
and /3 = 0.06. In Fig. 5, we show the probability for 
{foi fzOi fzzo} after fitting the model with the generated 
data. We have assumed three different hypotheses: the 
true values of {ot,j3}, {a = l,/3 = 0} and a “broad” 
marginalisation (a ~ _/V(l,0.05) and /? ~ iV(0.07, 0.05)). 
As shown in Fig. 5, the probability of fo is highly depen¬ 
dent on the choice of {a, /3} which may even lead to ruling 
out the true values of fo , when they are not known in ad¬ 
vance - as is the case when dealing with real data. The 
differences for the values of f z 0 and f zz 0 are less promi¬ 
nent. Note that the values for the “broad” marginalisa¬ 
tion do not cover the true values. However, the errors 
are very large for every case, which leads to a completely 
degenerated fit, such that a wide range of completely 
different f(R) models lie in the ler region. For exam¬ 
ple, by considering the viable Hu-Sawicki model [15], it 
is straightforward to show that the {fo, f z 0 , fzzo} of this 
model, for a wide range of random values of its free pa¬ 
rameters, lie within the ler region in Fig. 5, despite the 
fact that the mock data was generated by a very different 
model. Consequently cosmography is not able to distin¬ 
guish between them. Hence, if one decides to marginalise 
considering all possible values for a and /3, no useful con¬ 
straints can be derived. 


On the other hand, priors can bias the results strongly, 
which may be thought of as an avoidable limitation for 
current data given the size of errors, but may become 
an issue when hundreds of thousands of supernovae have 
been observed. Therefore, we see that the cosmographic 
approach is extremely weak at reconstructing f(R) grav¬ 
ity theories and therefore unable to rule out radically dif¬ 
ferent kinds of f(R) theories. The same issue arises when 
dealing with theories which have higher-order derivatives 
in the matter sector. In order to show this, let us just 
obtain the expansions within the simplest Galileon model 
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Figure 3: Left panel: Posterior probability for jo considering 4 parameters (0i), 5 parameters ( 6 1 ) and an XCDM model. Right 
panel: Posterior probability for the dark energy equation of state parameter w. 







Figure 4: Left upper panel: Scalar potential V{(j> ) evolution in terms of the redshift for four (#i) parameters expansion. The 
light/dark blue regions represent the la/2a regions. Right upper panel: Hubble parameter for model (4) as obtained from 
the fourth-order cosmographic parameters obtained in the MCMC analysis. Again light/dark blue regions represent the 1 <t/ 20- 
regions. Lower left and right panels are analogous to the upper counterparts for five (# 2 ) parameters expansion. 
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Figure 5: Effects of the marginalization for different choices of the free parameters a and /3 for f(R) theories. 


[16]: 


5 = 


d 4 xy/=g 


R + ~ (ci<j) + c 2 d fl ct)d ll (j) 


+ c^d^d^cjGcj)) 


(9) 


where Ci— 1 . 2,3 are coupling constants. The field can be 
expanded as 


Cl 

6 Hi 

4>zz0 

2</>zo 


4>0 = — 1 — «G — Pg + 


1 

"fc 


o a G o 

9 o - 3--2 

PG 


1 + -(1 + 90) 


( 10 ) 


where qg = C24>o z /6 and /3g = C 3 .R 0 </>[],. Thus, the 
same issues as in f(R) gravity appear, due to two ad¬ 
ditional free parameters ac and /3g- The only advan¬ 
tage Galileons theories have relate to the errors which 
are not as large as in the f(R) gravity case, since the 
gravitational sector does not include higher-order deriva¬ 
tives. Consequently the expansion of the scalar field up 
to second order only depends on go- In conclusion, the 
reconstruction process based on cosmography for theo¬ 
ries containing extra degrees of freedom - not only f(R) 
gravities or Galileons but also other extended theories as 
for instance, Horndeski-like Lagrangians - leads to com¬ 
pletely unconstrained parameters of the models under 
consideration. 


V. CONCLUSIONS 

In this paper, we investigated several issues that 
emerge when using the cosmographic approach to con¬ 
strain or reconstruct cosmological models, which have 
been overlooked in the previous literature. We first com¬ 
pared the cosmographic approach using different inde¬ 
pendent variables and found that when using the so- 
called y-redshift, biased constrains and larger errors 
are found on the cosmographic parameters compared to 
usual redshift z, even though y appears to be theoreti¬ 
cally better motivated. 

We then looked at what would happen if we used a 
fiducial model slightly different from ACDM, expecting 


that the fit of the cosmographic parameters would show a 
clear deviation away from the jo value in ACDM (jo = 1). 
Instead we found that the results depend considerably 
on the order of the cosmographic expansion, leading to 
ambiguous constraints (as shown in Fig. 3). 

Finally, and motivated by previous works in the liter¬ 
ature, we attempted to use cosmographic methods to re¬ 
construct different models of dark energy. We found that 
when considering theories with no higher-order deriva¬ 
tives in the gravitational or matter sectors, the method 
gives a clear picture of the underlying model, even though 
the errors are usually larger than those obtained using 
other approaches. However, when dealing with theories 
whose Lagrangians contain higher-order derivatives in ei¬ 
ther the gravitational or matter sector, there are extra 
free parameters which cannot be constrained by cosmog¬ 
raphy and must be marginalised in order to recover the 
underlying theory. As we illustrated, this fact leads in 
general to large errors preventing us from ruling models 
out with very different cosmological background evolu¬ 
tion. Therefore, given the current state of affairs, it seems 
that the cosmographic approach is not a useful tool for 
theory reconstruction. 

The focus here has been to highlight a number of short¬ 
comings of the cosmographic approach, but we should 
also ask if there are other limitations that need to be 
taken into account and whether it is possible to overcome 
the limitations discussed in this paper. Regarding the 
first question - there are two effects which were not con¬ 
sidered in our analysis. The first one relates to the role of 
spatial curvature in the constraints, since it is known that 
Hfc ^ 0 can induce a time variation for the dark energy 
equation of state [17]. Also, effects such as gravitational 
and Doppler lensing [18], or even local gravitational red- 
shifts [19] lead to an extra scatter in the Hubble diagram, 
which degrades cosmological constraints when using SNe 
la data. Although quantified in a ACDM model using 
simulations and inserted as an extra error in SNe la anal¬ 
ysis [9] , the impact on extended gravity theories may be 
more important and a thorough study of this effect is 
necessary in order to rule out those theories for which 
the cosmographic approach cannot be applied. 

Finally in response to the second question - we can 
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identify a number of challenges faced by this method: 
There needs to be i) a clear definition of auxiliary vari¬ 
ables, together with their range of applicability and ex¬ 
tensive testing against mock data; ii) a robust statistical 
method which could in principle establish a trade-off be¬ 
tween number of data points, number of cosmographical 
parameters and Bayesian evidence in order to provide 
criteria to safely rule out models; and in) very well mo¬ 
tivated priors over the extra parameters in higher-order 
theories in order to get useful constraints for these theo¬ 
ries, even when neglecting lensing and local effects. Only 
when these issues are addressed, can we be confident that 
cosmography is a useful tool for constraining extended 
theories of gravity. 
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